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We consider the Bianchi identity as a field equation for the distortion of the shapes of images 
produced by weak gravitational lensing. Using the spin coefficient formalism of Newman & Penrose 
Q, we show that certain complex components of the Weyl and Ricci curvature tensors are directly 
related to fundamental observables in weak gravitational lensing. In the case of weak gravitational 
fields, we then show that the Bianchi identity provides a field equation for the Ricci tensor assuming 
a known Weyl tensor. From the Bianchi identity, we derive the integral equation for weak lensing 
presented by Miralda-Escude 0, thus making the Bianchi identity a first principles equation of weak 
gravitational lensing. This equation is integrated in the important case of an axially symmetric lens 
and explicitly demonstrated in the case of a point lens and a SIS model. 
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^ ■ I. INTRODUCTION 

The Bianchi identity in general relativity is usually expressed as 



> 
00 



V [o i? H / = (1) 



and is considered to be a consequence of the geometrical properties of the covariant derivative operator. Most 
students of relativity come across the Bianchi identity in a contracted form that is quite useful in discussions of 
\^ , various properties of the Einstein tensor. 

' The decomposition of the Riemann tensor into the Ricci (R a b) an d Weyl (C a b C[ j) tensors, 

m : 

O . Rabcd = Cabcd + (g a [ c Rd]b - 9b[cRd]a) - -^Rg a [c9d]b, (2) 

cr 6 

, presents a powerful use of the Bianchi identity. As Hawking & Ellis 3] point out, the Ricci tensor is determined by 
the Einstein field equation, while the Weyl tensor is the part of the curvature not determined locally by the matter 
distribution. However, the application of the Bianchi identity to Eq. [21 yields a constraint on the Weyl tensor once 
• *h . the Ricci tensor is found from the Einstein field equations. By treating the Bianchi identity as a field equation for the 
Weyl tensor given a known Ricci tensor, Newman & Penrose pj (hereafter NP) developed a consistent methodology 
^ , for the study of gravitational radiation. 

In this paper, we consider the Bianchi identity in the reverse manner. Specifically, we assume that the Weyl tensor 
is known and treat the Bianchi identity as a partial differential equation for the Ricci tensor. 

Our motivation for this approach is the study of weak gravitational lensing, where measurements are made of 
the distortion of images due the presence of gravitational fields. From the observed distortion of images, a matter 
distribution can be inferred, so weak gravitational lensing is now a very active area of observational research with 
many papers published each year. Pivotal early work relating the observed distortion and inferred mass distribution 
was done by Miralda-Escude and Kaiser & Squires Q amongst others. Substantial reviews of the subject are 
provided by Mellier |(| and Schneider & Bartelmann 0. For this paper, we are most interested in the presentation 
given in Miralda-Escude [2j. 

In this paper, we present four results. First, we show that the observed distortion and inferred matter distributions 
of weak lensing are equal to projected components of the Ricci and Weyl tensors as expressed in the formalism 
introduced by Newman & Penrose Q]. Next, we demonstrate that the integral equation used in observational weak 
lensing to determine the projected mass distribution from observed image distortions is actually an integral relation 
between the >Fo and <&oo components of the Weyl and Ricci tensors. 
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Our most important result is that we prove that a particular component of the Bianchi identity yields a partial 
differential equation for the $00 component of the Ricci tensor given a known 'Fq Weyl tensor. Using a Green's 
function, we derive the standard equation used in observational weak lensing as formulated in Miralda-Escude Q 
from the Bianchi identity. Thus, we have developed a PDE approach to weak lensing that may provide a new 
calculational approach to the field that allows for more accurate determinations of mass mappings. 

Finally, we consider the important case of axially symmetric lenses and show that in this case the Bianchi identity 
is easily integrated. Two explicit examples are shown: a point lens and a singular isothermal sphere. 

Our results on weak lensing extend a new way of thinking about gravitational lensing based on a space-time 
perspective. Strong lensing and general image distortion have been studied by several researchers including Frittelli 
& Newman || and two connected papers by Frittelli et al. 0,^3 with a substantial review provided by Perlick [Tl) . 
This paper is the first paper that seeks a first principles equation for observational weak lensing. 



II. COMPUTATION OF RICCI AND WEYL TENSORS 



In this section, we compute the complex Ricci and Weyl tensor components and the spin coefficients introduced in 
Newman & Penrose pj. Our goal is to determine the spin coefficient form of the Ricci and Weyl tensor components 
to first order in a static metric perturbation off flat space. 

While modern evidence strongly indicates that the universe is accelerating and dominated by dark energy, it is 
reasonable to chose to perturb off flat space for several reasons. First, we wish to directly compare our derivation of 
the equations of weak lensing with the standard derivations using the thin- lens approximation. Because all the "lensing 
activity" in the thin-lens approximation happens when light rays passes the lens, the choice of cosmology does not 
directly enter the equations. In fact, both in the literature and in actual application, the only place the cosmological 
model enters weak lensing calculations is in the conversion between redshift and angular diameter distance. Hence, one 
of the strengths of lensing studies is that the basic equations of the theory are relatively independent of cosmological 
model. 

Furthermore, since the basic cosmological models are conformally flat, our basic results can be translated to them. 
Beginning with a more "realistic" metric does not add to the calculation and introduces a number of conceptual issues 
best handled in a different paper. 

The physical situation that we wish to study is the appearance and distortion of a small patch of the sky containing 
extended galaxies that has been weakly lensed by some matter distribution. Both the astronomer in question and the 
sky of galaxies are assumed to be situated in reasonably isolated regions of space-time where the space-time is flat. 

The observed galaxies are connected to our astronomer by light rays along her past light cone. In fact, we may 
think of pencils of light rays that travel from each extended galaxy to the telescope. The path of these light rays is 
to be determined by integrating the null geodesic equations of the space-time metric, so that we may find a tangent 
vector, £ a , to the light rays. 

The matter distribution that lenses the light rays is encoded into a space-time metric. For our purposes, we may 
assume a weakly perturbed Minkowski metric, 

ds 2 = (1 + 2cp)dt 2 - (1 - 2tp){dx 2 + dy 2 + dz 2 ), (3) 
where ip is a static gravitational potential satisfying 



VV = 4np(x,y,z). (4) 

As the astronomer and observed galaxies are to be isolated, we assume that the matter distribution, p(x,y,z), is 
contained in some region of space far from the astronomer and galaxies. In practice, this is always the case. 
Following the program of NP, we chose a tetrad of null vectors, 



\? = (\ a u \ a 2 ,\ a 3 ,\Z) = (r,n a ,m a ,m a ), (5) 

associated with the pencil of light rays connecting our astronomer and each individual observed galaxy, where £ a is 
the real vector tangent to the pencil. By convention, the other vectors are chosen such that 



£ a n a = 1 & m a m a = -1, 



(6) 
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with all other products zero. m a and m° are spatial, complex null vectors parallel propagated along the pencil. 

In principle, the space-time metric can be written in terms of a general null tetrad as g ab — rf^XfXj, although this 
is not central to our presentation. Further, we will make physical restrictions on the tetrad vectors tied to our goal of 
finding the Weyl and Ricci tensor components to first order in the perturbation. 

We orient our spatial coordinates, (x,y,z), such that the astronomer is located at x — y = and z = z a with 
the "center" of the matter distribution at the spatial origin. With this placement, our astronomer's telescope points 
straight down the z axis. 

In principle, each null tetrad associated with an individual pencil of light connecting the astronomer to a lensed 
galaxy varies along the pencil. However, in our calculations, we are justified in considering each tetrad to be a constant 
null tetrad given by 



£" = -±=(1,0,0,1), n Q = -i= (1,0,0,-1) & m a = -i=(0,M,0), (7) 

for two reasons. First, the wide field telescopes used today see a very small portion of the sky. For example, the 
proposed 8.4 m Dark Matter Telescope would have a field of view of only 260 square milli degrees Hence, with 
our orientation, all the pencils of light are parallel to the z axis. Second, even though the individual pencils will be 
deflected by the lens, the deflection of the light ray will be proportional to the 2-dimensional (x, y) gradient of <f. 
Since we will be contracting the tetrad vectors with the Ricci and Weyl tensors, to work to first order in y>, we must 
neglect any variation in the null tetrad. 

The basis of the NP formalism is twelve complex spin coefficients that are analogous to the twenty-four real Ricci 
rotation coefficients defined by 

Yjk = A"A fc VbA^. (8) 

Because our choice of tetrad involves vectors whose components are all constant, all the NP spin coefficients are zero 
for our physical situation (to zeroth order in the metric perturbation). 

In the NP formalism, tetrad components of the Ricci and Weyl tensors are computed by contracting the coordinate 
components with the tetrad vectors. This yields five complex Weyl tensor and ten complex Ricci tensor components 
(see appendix 0. 

Using the null tetrad above, the non-zero NP formalism Ricci curvature components are 

$oo = $22 = 2*n - ^VV (9) 

The first order Weyl tensor components are 

*0 = *4 = ^ (fxx - Vyy + 2i<Pxy) > 
*l = *3 = ~ {fxz + i^yz) , 

*2 = \ Uzz - gVV) • (10) 

We note that up to numerical factors, <E>oo is equal to the matter density p(x, y, z) that determines the gravitational 
perturbation ip. 

The four directional derivatives associated with the null tetrad play an important role in the NP formalism and are 
given special names: 

D = £ a d a , A = n a d a , 6 = m a d a , k 6 = m a d a . (11) 
Using our choice of tetrad, and the fact that the metric perturbation is static, these derivative operators will act as 



D 



i e_ 

V2 dz> 



5 ~ V2 \dx +l dy) ' 



A 



6 = 



1 d_ 

V2 8*' 



J_ (JL- sJL 

y/2 \ox dy 



(12) 
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III. OBSERVATIONAL WEAK LENSING IN PRACTICE 



In this section, we briefly outline the usual presentation of the equations of weak lensing by a thin lens so that 
we can draw parallels to our presentation that uses the Weyl and Ricci tensor. The derivation presented here relies 
heavily on that of Miralda-Escude Q • 

The practicing astrophysicist assumes a thin (or two dimensional) lens that divides a background space-time into 
an observer and a source side. The lens lies in a two (spatial) dimensional "lens plane" that is perpendicular to the 
line of sight of the telescope. The observed background galaxy is said to lie in a two dimensional "source plane," also 
perpendicular to the line of sight. 

Using dimcnsionless cartesian coordinates f a in the source plane and f\ = (x,y) for the lens plane, the starting 
point for thin-lens lensing is a mapping from the lens plane to the source plane given by 

r s =n-^2Wi)- (13) 

The two-dimensional gravitational potential tp is determined from an ordinary three-dimensional potential ip by 
projection, 



ip(f f i)=J dzcp(x,y,z), (14) 

where z a and z s are the z coordinates of the observing astronomer and distant source, respectively. The two dimen- 
sional potential ip will satisfy 



1^ = ^ + ^=^, (15) 

where £(rj) is a projected mass density and S cr i t is a critical density for strong lensing, or the appearance of multiple 
images. 

To consider weak lensing, one forms a Jacobian matrix 



dr s _ ( 1 - ip xx -ip xy \_{1-k-\ -p 



J =-zr- = [ x ^ XX i 7 = i 1 x > ( 16 ) 

n V -1p xy 1 - Ipyy J V ~ fJ, 1 - K + A ' 



where 



\ {i>XX + *Pyy) = ^A 2 ^ = (17) 



is referred to as the "convergence," and the two quantities 



A = - (l/j xx - Ipyy) & M = i>xy (18) 

are called the "shears" in the thin-lens literature. Frittelli et al. 01 shows how these quantities are related to the 
convergence (p) and shear (a) in general relativity. For this paper, it is critical to note the similarities that k has 
with $00 and that ^o has with A and /i. 

If the outer surface of a circular, extended source is parameterized by t s = R (cos 8, sin 5), then under the thin lens 
mapping the inverse of the Jacobian matrix in Ea. 1161 maps t s —> ti, a new elliptical curve in the lens plane. The 
orientation and ellipticity of the resulting ellipse is determined by A and p, - see Schneider et al. 0] for the details. 

Since the orientation and ellipticity is observable, the goal in observational weak lensing is to measure the shears 
(A, p) and infer the projected mass density, K. To do this, one inverts Eq. El using the two dimensional Green's 
function and applies the differentiations in Eq. ^|to the result. This yields 



A 



dr 



K>(f) (— cos2?7) 



7T 



K (^) (— sin2?7) 



(19) 
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where rj is the angle in the lens plane between f '— f 1 and the x axis. 

By employing a Fourier transform technique, Miralda-Escude shows that one can invert Eq. ^|to obtain 

« ( f) = / # m>m . [£g^M, (20) 



Equation 1201 is the primary equation for observational weak lensing. It is an integral relation that allows one to 
determine the projected mass density at every point if one has measured the "shears." 

IV. WEAK LENSING WITH CURVATURE TENSORS 

We may exploit the similarities between the definitions of the thin-lens convergence and shears and <E>oo and to 
recast Eq. [2U]in different language. The first step is to project our gravitational potential and curvature tensors into 
a "thin-lens" form. Equation 1141 projects our metric perturbation tp(x,y,z) to a two dimensional ip(x,y). Then we 
define 

1 

L $00 = J dz $00 = ^ {ipxx + ^yy) , (21) 

where we use Eq. 03 with Eg. 1141 and the property that tp z is zero far from the lens. Likewise, we can define 



L *0 = y dZ^Q= -{l(> XX -i>yy + 2ilp x y). (22) 

We note that l$oo = k and = A + ifi, so that the fundamental observable in weak gravitational lensing is 
a projected component of the Weyl tensor, while the fundamental inferred quantity of observational interest (the 
projected mass density) is a projected component of the Ricci tensor. 

Then exactly as in section [TTTl one can use a two-dimensional Green's function to invert Eq.[2]for ip and plug this 
relation into Ea. l22l to obtain an integral equation specifying given a known l$oo- Taking the Fourier transform 
yields a relation between the transforms of an d l&oo which is easily rearranged such that the inverse Fourier 
transform produces the desired integral equation for l&oo, 



L $oo (r) = - / dr 1 ™li , (23) 

J tt \r — r\ z 

where 77 again is the angle between r — r* and the x axis in the lens plane. 

Equation[2Bldemonstrates that the fundamental equation of weak lensing is an integral relation between components 
of the Weyl and Ricci curvature tensors. Specifically, a projected version of $0 is known through observation and 
used to determine a projected mass density l$oo- 



V. THE BIANCHI IDENTITY 



As discussed in section applying the Bianchi identity, Eq. ^ to the decomposition of the Riemann tensor into 
the Ricci and Weyl curvature tensors produces a differential constraint on the Weyl tensor. To obtain the Bianchi 
identity in the NP spin coefficient formalism, one contracts the Bianchi identity with all possible combinations of 
the tetrad vectors and writes the resulting equations out using the definitions of the spin coefficients and Ricci and 
Weyl tensor components. In the non-vacuum case, this results in twelve equations, which are listed in a number of 
references including Newman & Tod |l4j |. 

The first full Bianchi identity equation is 



- D*i + L»$oi - <S$oo = (4a - 7r)* - 2(2p + e)*i + 3k* 2 + {n-2a- 2/3)<J> o 

+ 2(e + p)$ i + 2o-$io - 2k$h - k$ 2- (24) 
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Due to our choice of tetrad, all the spin coefficients are zero, so the entire right hand side of Eg. 1241 vanishes. Also, 
$01 is zero, so to first order in the gravitational perturbation, we have 

<J* - D^i - 5$ 00 = 0. (25) 

Equation 1251 is a partial differential equation that holds at every point along the pencil of rays connecting our 
astronomer with the distant observed galaxy. To compare this relation with Eq. [23J we project into a lens plane by 
integrating out the z coordinate. Integrating in z from z a to z s effectively yields zero since D ~ d z and "fi 

evaluated far from the lens is assumed to be zero. Passing the z integration through the 5 derivative operators, we 
have 

6l^oo^S l ^o- (26) 

We consider Eq.[2H|as a field equation for l^oo given a known l^o- Equationl26lrelates the fundamental quantities 
of weak gravitational lensing and is derived from first principles in general relativity. 

VI. GREEN'S FUNCTION 

In this section, we show that the relation between the Weyl and Ricci tensor derived from the Bianchi identity, 
Eq. [23 is the differential version of the integral relation used in standard weak lensing studies, Eq. [231 To do this, 
we employ a set of Green's functions first developed by Porter 0] and expanded in Ivancovich et al. which we 
refer to as the Porter Green's functions. 

The Porter Green's functions are designed to work with powers of the S and 9 differential operators. The S derivative 
operator acts on functions on the sphere of different spin weight in different ways and serves as a spin- weight raising 
and lowering operator Functions on the sphere of a spin weight s are expandable in a series of spin- weighted 

spherical harmonics s Yi m . Because different powers of 3 and § annihilate the s Y/ m with \s\ = I, a set of Green's 
functions can be developed. 

Our main result, Eq. 1261 is a partial differential equation that holds in an (x, y) plane that we consider to be the 
lens plane. To make contact with the Porter Green's functions, we introduce the pair of complex coordinates (£, Q 
defined by 

( = -L(x-iy). (27) 

In these coordinates, m a points in the £ direction and 8 — 8q. 

By stereographic projection of the sphere into a equatorial plane from the pole, the (£, C) coordinates are complex 
coordinates for the sphere. With the exception of the point at infinity, where the projected curvature tensors 
and l$oo are assumed to be zero, introduction of the complex stereographic coordinates turns Eq. 1261 into a partial 
differential equation on the sphere. 

Since l$oo is a spin- weight zero function, the application of S to it will take the form 

3l*oo = + (28) 
For this reason, we multiply both sides of Eq. 1261 by (1 + ££) and our Bianchi identity takes the form 

8 £ *oo = (l + CC)*£*o=Ai(C ) C), (29) 

where Q is a function of spin-weight 1. 

The Porter Green's function for an equation of the form in Eq. [551 is 

l$oo(C,C)= / i^ ,-i(C,C;^»7)(l + w)^i*orf|«n' ( 3 °) 
Js 2 

where the integral is taken over the sphere parameterized by complex coordinates (r], fj) with the area element 
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2 dr\ A dr\ 

djJLr, - 



i (1 + r]r]) 2 

The kernel of the Green's function for an equation of the type in Eq. is 



K 0t _ 1 ((,{;r 1 ,fi) = ±l±^>. (31) 
4?r C — V 



Putting all this together, we have 

i_ r s vL Mv) 

Hit J S 2 (-fj 

It is convenient to multiply by 1 in the form 1 = (C ~ v)/(C ~ 7 l) to obtain 



l$oo(C) = ^/ 7 vy _" drjAdfj. (32) 

2lTT J S 2 Q-7] 



l$oo(C) = ^r- I C * rS v L*o(v)dv A dfj. (33) 

2«7i- J s * (C-v)(C-v) 

Because we want to show that Eq. 123 is equivalent to Eq. [231 we reintroduce cartesian coordinates 

( = -j=(x-iy) rj = -j=(x'-iy'). (34) 

In these coordinates, Ea. 1331 becomes 

l**>{x, y) = ^J [ £Z*)2~+( y Z$ - **) d *' d y'- ( 35 ) 

Comparing Eq. 1351 with Eq. 1231 we see that we need to flip the 5^ derivative operator that is inside the integral by 

integrating by parts. Using the definition of the angle r\ as the angle the vector r — r 1 makes with the +x axis and 
integrating Eq. 123 by parts yields 

l$oo(^,2/) = - / dr — — —2 (36) 

J 7r \r — r\ A 

plus the surface term 

I J <„w ^ - «v) { ifezf^pia . (3r , 

One can see that the surface term will vanish by breaking the integral into two parts for the d x i and d y i derivatives. 
Each part vanishes as the integrand is evaluated at either x' — oo or y' = oo. 

Therefore, we have shown that the application of the Porter Green's function to the Bianchi identity, Eg. 1261 results 
in the integral relation used in weak gravitational lensing studies. 



VII. AXIAL SYMMETRY 



In this section, we consider those mass distributions which are axially symmetric around the line of sight. In practice, 
this is a very important simplifying assumption used in gravitational lensing studies. Under this assumption, Eq. 1261 
is easily integrated directly, without using the Porter Green's function. 

First, we change coordinates in the lens plane from the cartesian (x , y) coordinates to standard axial coordinates 
(r , 0). In the axial coordinate system, the 6 directional derivative operator in Eo. 1121 is 
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5 = 71 [ dr + l d V ■ (38) 

As an ansatz for an axially symmetric lens, we assume that the Ricci tensor is a function of r only, 

l$oo= l^ooW, (39) 

and that the Weyl tensor can be written as 

L * = T(r)e 2 ^. (40) 

This second assumption is motivated by the concept of spin weight in the NP formalism (l4| . 

Applying Eq. |2Hlto our ansatz for the solution to Eq.l^Hl the </> dependence drops out, and one is left with 

d rL ^ 00 (r)=d r T(r) + -T(r), (41) 
r 

which is now a one dimensional equation. 

Recall that we have shown that l^o or Y(r) is a known, measurable quantity and that l$oo is essentially the 
projected matter density. This means that after integrating over r, Eq. 1411 gives us the functional form of the 
projected matter distribution: 

L $ooW = T(r) + J^-Tdr + C. (42) 

The constant of integration represents the mass sheet degeneracy present in gravitational lensing and can be taken as 
zero. 

VIII. TWO AXIALLY SYMMETRIC EXAMPLES 

In this section, we examine two axially symmetric mass distributions in the context of weak lensing as expressed 
by Eg. 1421 The models we consider include a singular, Schwarzschild-like point lens and a singular isothermal sphere 
(SIS) model. For each model, we will determine the gravitational potential, ^p{r) 1 then compute the projected Weyl 
and Ricci tensor components and show that they obey Eq. 14*2*1 

A. Point lens 

As our first example, we consider a Schwarzschild point lens of mass M at the origin. For this mass distribution, 
the gravitational potential i\) is given by 

tp = 2M ln(r) = 2M In (^/x 2 + y 2 ) . (43) 
One can show by direct computation that 



^ = _^e 2 ^, (44) 



which matches Eq . 1401 yielding 
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Inserting this functional form for T into the right hand side of Eq. 021 gives 

-2M f , -AM 

— — + / dr — — = 0. (46) 

Equation 1461 indicates, through the Bianchi identity Eg. 1421 that the projected Ricci tensor is zero at all radii r > 0, 
which is the expected result for a Schwarzschild, point-like lens. 



B. SIS model 



The SIS model is a one parameter model given by 



£(r) 



2r' 



(47) 



where a v is the velocity dispersion. Although unphysical due to the singularity at the origin and infinite total mass, 
SIS models do explain the flat rotation curves of galaxies and are widely used as models in gravitational lensing 
studies. 

For our purposes, the simplest way to integrate the projected gravitational potential for the SIS model is to follow 
the argument of Schneider et al. [13| in section 8.1, where it is shown that one can write the projected gravitational 
potential as 



ip(r)=8n I r' dr'S(r') In (-) . (48) 



It is then relatively simple to show that 



ip = A-iral r > (49) 

so that the projected Weyl tensor is 

L * = -^ e «f (50) 
r 

From Eq. 0OJ we have T = —2ira1/r and from the Bianchi identity 

r + [ dr ™ = 2 -^. (5i) 

J r r 

By direct computation, 



L$00 = \ {^xx + Ipyy) = (52) 



so that the Bianchi identity, Eq. 021 is preserved. 



IX. DISCUSSION 



The main results of this paper are that the fundamental quantities of weak gravitational lensing are projected 
components of the Ricci and Weyl tensors and that the Bianchi identity provides a field equation for the projected 
matter density derived from first principles. We explicitly show that the common integral equation used in weak 
gravitational is an integral version of the Bianchi identity. To our knowledge, these results are known in the lensing 
or relativity communities. 
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From this perspective, this paper extends the work presented in a series of recent papers that has helped reunite 
applied gravitational lcnsing with its roots in general relativity. The underlying perspective of these papers has been 
that any observed lensing phenomena must be encoded into the past light cone of the observer. 

A primary purpose of these papers has been to identify the equation from general relativity that is approximated in 
the thin lens treatment used by practicing astrophysicists. Frittelli & Newman 8] first pointed out that the lens should 
be coded into a space-time metric, and then solving the null geodesic equations would simultaneously give the "time 
of flight" equation and lens mapping. In a two-paper set, Frittelli et al. Q related general relativity's optical scalars 
to the "shears" and "convergences" cited in the thin-lens literature. These papers also wrote out integral relations 
that showed how image distortion grew continuously along the pencil of rays connected a source and observer. A 
generalization of the Fermat principle of least time was present in Frittelli et al. ^3 • 

Because Eg. 1231 is kincmatically derived, it does not indicate a true connection between observational weak lensing 
and relativistic first principles. The derivation of Eq. 1231 relies only on the definitions of $oo an d ^o- 

However, we show in this paper that the projected Bianchi identity, Eq. is t ne differential version of Eq. [23 
Thus, we have a first principles derivation of the basic equation of weak gravitational lensing for the first time and a 
potentially useful new PDE approach to the subject. 

Most directly, Eq. 1261 provides a new PDE approach for determining the projected mass density in weak lensing 
studies. While numerical approaches to solving PDEs with real data tend to less stable than integral approaches, they 
also tend to provide higher accuracy, so it is not unreasonable to pursue approaches to data analysis based directly 
on these results. 

Another potential application of this work might be in the cumulative analysis of thick lenses or repeated lensing 
by multiple clusters along the line of sight. In this case, one would not apply a thin-lens assumption, but would use 
Eq. [23 as a field equation along the line of sight for a full metric that encodes a more realistic cosmological model. 
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APPENDIX A: WEAK FIELD RICCI AND WEYL TENSORS 

In this appendix, we give the usual component form of the Ricci and Weyl tensors considered in the text, along with 
their definitions in the Newman & Penrose [l( formalism. In Eqs. IA1IA3I the usual Einstein summation conventions 
are not employed. 

For the weak field metric in Eq. |3| to first order in the perturbation, tp, the non-zero Christofel symbols are 

r° 0l = & r* ife = -tp k 

r\x> = w r i fcfc = Vi (i^fe) (Ai) 

where tp+ = di<p, denotes time and i, j, k denote spatial components. The non-zero, first order Ricci tensor components 
are 

R 00 = -VV R« = -VV (A2) 
with V 2 = <9 2 + dy + 9 2 , while the non-zero, first order Weyl tensor components are 

Com = ^ (-3<£>ii + V 2 </j) Coioj = -<fij i ^ j 

C m = 3 (3<Pkk - V 2 <£>) i ^ j ^ k Cijik = -ifjk j ^ k. (A3) 

In the NP formalism, the components of the Ricci and Weyl tensors are contracted with the null tetrad to give 
named components. The Weyl tensor components are defined by 

*o = -C abcd £ a m b £ c m d , *i = -C abcd £ a n b £ c m d , 
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*2 = -~ {C abcd tn b l c n d - C abcd tn b m c m d ) , 

$3 = C abcd tn b n c rh d , $4 = -C abcd n a fh h n c m d , (A4) 
and the Ricci tensor components are 

$00 = -\R ab e a e b , $10 = -\R ab l a fh b , $ 20 = -^R ab fh a m\ 

$01 = -^R ab i a m\ $n = -\ (R ab e a n b + R ab m a m b ) , $ 2i = ~R ab n a m b , 

$02 = ~R ab m a m b , $12 = -\R ab n a m b , $ 22 = -^R ab n a n b , 

A = ^i?. (A5) 
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